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We demonstrated theoretically that the renormalization of the electron energy spectrum near 
the Dirac point of graphene by a strong high-frequency electromagnetic field (dressing field) dras¬ 
tically depends on polarization of the field. Namely, linear polarization results in an anisotropic 
gapless energy spectrum, whereas circular polarization leads to an isotropic gapped one. As a con¬ 
sequence, the stationary (dc) electronic transport in graphene strongly depends on parameters of 
the dressing field: A circularly polarized field monotonically decreases the isotropic conductivity of 
graphene, whereas a linearly polarized one results in both giant anisotropy of conductivity (which 
can reach thousands of percents) and the oscillating behavior of the conductivity as a function of 
the field intensity. Since the predicted phenomena can be observed in a graphene layer irradiated 
by a monochromatic electromagnetic wave, the elaborated theory opens a substantially new way to 
control electronic properties of graphene with light. 


INTRODUCTION 


Since the discovery of graphene [l| , it has attracted the persistent interest of the scientific community. Particula^, 
the influence of an electromagnetic field on the electronic properties of graphene is in the focus of attention [2, Q . 
Usually, the electron-field interaction is considered within the regime of weak light-matter coupling, where the electron 
energy spectrum is assumed to be unperturbed by photons. However, a lot of interesting physical effects can be 
expected within the regime of strong light-matter coupling, where the electron energy spectrum is strongly modified 
by a high-frequency electromagnetic field. Following the conventional classification, this regime is jurisdictional to 
quantum optics which is an established part of modern physics ii- Therefore, the developing of interdisciplinary 
research at the border between graphene physics and quantum optics is on the scientific agenda. 

The methodology of quantum optics lies at the basis of various exciting fields of modern physics, including quantum 
information Q, polaritonics Q, quantum teleportation ii, quantum cryptography [a, El, etc. Particularly, 
it allows to describe fundamental physical effects (e.g., Bose-Einstein condensation of polariton s Il2ll and optical 
bistability 131) and creates a basis of modern technological applications (e.g., optical logic circuits [1^, novel sources 
of terahertz emission 151, and novel types of lasers [1^ llll)- Within the quantum optics approach, the system 
“electron -|- strong electromagnetic field” should be considered as a whole. Such a bound electron-field system, which 
was called “electron dressed by field” (dressed electron), became a commonly used model in modern physics ii- 
The field-induced modification of the energy spectrum and wave functions of dressed electrons was discovered many 
years ago and has been studied in detail in various atomic systems [isl - Esjl and condensed matter 24 - ^ . In graphene- 
related research, the attention has been paid to the field-induced modification of energy spectrum of dressed electrons 
341441|. optical response of dressed electrons [42|, transport of dressed electrons in graphene-based p-n junctions [43 1 
and electronic transport through dressed edge states in graphene 44-4^. As to stationary (dc) transport properties 
of a spatially homogeneous graphene layer dressed by light, they still await detailed analysis. The present Report is 
aimed to fill partially this gap at the border between graphene physics and quantum optics. 


MODEL 


For definiteness, we will restrict our consideration to the case of electron states near the Dirac point of a single 
graphene sheet subjected to an electromagnetic wave propagating perpendicularly to the graphene plane. Let the 
graphene sheet lie in the plane {x, y) a.t z = 0, and the wave promgate along the z axis [see Fig. (|T|)]. Then electronic 
properties of the graphene are described by the Hamiltonian [2,0 

% = vcr ■ (fik — eA), (1) 

where cr = (era,, Gy) is the Pauli matrix vector, k = {kx, ky) is the electron wave vector in the graphene plane, v is the 
electron velocity in graphene near the Dirac point, e is the electron charge, and A = {Ax, Ay) is the vector potential 
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FIG. 1: Sketch of the electron-field system under consideration: The graphene sheet dressed by (a) circularly polarized elec¬ 
tromagnetic wave with the amplitude Eq and (b) linearly polarized one. 


of the electromagnetic wave in the graphene plane. In what follows, we will be to assume that the wave frequency, 
w, lies far from the resonant frequencies of graphene, 2vk. Solving the non-stationary Schrodinger equation with the 
Hamiltonian 

= Uipk, 
ot 

we can obtain both the energy spectrum of electrons dressed by the electromagnetic field, Ek, and their wave functions 
'0k (see technical details of the solving within the Supplementary Information attached to the Report). 

For the case of the circularly polarized electromagnetic field with the vector potential 

/ ifo . ^0 ■ , 

A = — cos U!t, — sinwr 

V W w 

we arrive at the energy spectrum of the dressed electrons, 

Ek = ±sJ{£g/2Y + {hvkY, (2) 

where signs “-I-” and ” correspond to the conduction band and valence band of graphene, respectively. 

Eg = ■sj {hu))^ + {2veEo/u!)'^ — hu! (3) 

is the field-induced band gap in graphene, Eq is the amplitude of electric field of the electromagnetic wave, and the 
field frequency cu is assumed to satisfy the condition of a; ^ ^2veEo/h. Corresponding wave functions of electrons 
dressed by the circularly-polarized field read as 


0k = V5k(r)e 


—ie-k.tjh 


Ek|TEg/2 ^_,g/2 ^ / kk|±Eg/2 ^,g/,^^^^^ 


2kk| 


2kk| 


(4) 


where r = {x,y) is the electron radius-vector in the graphene plane, $i_ 2 (r) are the known basic functions of the 
graphene Hamiltonian (the periodical functions arisen from atomic 7r-orbitals of the two crystal sublattices of graphene) 
y, i^k(r) = /\fS is the plane electron wave, S is the graphene area, and 9 is the azimuth angle of electron in the 
space of wave vector, k = (fccos0, fcsinS). 

In the case of linearly polarized electromagnetic field with the vector potential 

A = coswf, 0 

directed along the x axis, the energy spectrum of the dressed electrons reads as 

Ek = ±hvkf{9), (5) 

and the corresponding wave functions of dressed electrons are 


sin0 


/ 2veEo \ 


0k = I [<hi(r) ± $2(0] X - 0osg + /(g) ^° ^ 


^^i{veEo jhu)^) sin cut 




( 6 ) 
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where 

f{9) = ^ cos^ 0 + Jo(^ ° ^ sin^, (7) 

Jo{z) is the Bessel function of the first kind, and the field frequency to is assumed to satisfy the condition of a; ^ vk. 



FIG. 2: The energy spectrum of dressed electrons in graphene for the dressing field with different polarizations: (a) circularly 
polarized dressing field; (b) dressing field polarized along the x axis. The energy spectrum of electrons in absence of the dressing 
field is plotted by the dotted lines and ef is the Fermi energy. 

The energy spectra of dressed electrons, m and ®, are pictured schematically in Fig. [2l As to a consistent 
derivation of Eqs. @-0, it can be found within the Supplementary Information attached to the Report. In order 
to verify the derived expressions, it should be stressed that the energy spectrum of electrons dressed by a classical 
circularly polarized field, which is given by Eq. exactly coincides with the energy spectrum of electrons dressed by 


a quantized field in the limit of large photon occupation numbers 36| . This can serve as a proof of physical correctness 


of the presented approach elaborated for a classical dressing field. 

In order to calculate transport properties of dressed electrons, we have to solve the scattering problem for nonsta¬ 
tionary electron states ([4]) and ([H). Eollowing the scattering theory for dressed conduction electrons [s^, the problem 
comes to substituting the wave functions of dressed electrons O and ([6]) into the conventional expression for the Born 
scattering probability [d^. Assuming a total scattering potential in a graphene sheet, U{r), to be smooth within an 
elementary crystal cell of graphene, we can write its matrix elements as 

($i(r)(/)k'(r) |C/(r)| $j(r)(/:k(r)) ~ t/k'k%, 

where C/k'k = (</^k'(r) |t^(r)| (y3k(r)), and Sij is the Kronecker delta. As a result, the Born scattering probability for 
dressed electronic states in graphene takes the form 


2vr Q 2 

U'k'k = ^IXk'kl |C7k'k| i5(£k' — Ek), 
h 


( 8 ) 


where 


/|ek'|-£g/2 /|ek|-eg/2 i{0'-e)i2 , /kk'|+^s/2 / kk| + eg/2 

2|£kH V 2|£k| " 2|ek/ V 2|ek| " 


(9) 


for the case of circularly polarized dressing field, and 


Xk'k = 


I cos{9') + f{6') cos{9) + f{9) 


2/(d') 


2/(0) 


sin(d') sin(0) 2 /^ 2 uei?o\ 

+ (00 + /(0O cos(0) + f{9) ) 


COS 


( 10 ) 


for the case of linearly polarized dressing field. 

In what follows, we will assume that the wave frequency, w, meets the condition 


WTo > I, 


( 11 ) 
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where tq is the electron relaxation time in an unirradiated graphene, which should be considered as a phenomenological 
parameter taken from experiments. It is well-known that the intraband (collisional) absorption of wave energy by 
conduction electrons is negligibly small under condition m (see, e.g., Refs. [H S [HI). Thus, the considered 
electromagnetic wave can be treated as a purely dressing field which can be neither absorbed nor emitted by conduction 
electrons. As a consequence, the field does not heat the electron gas and, correspondingly, the electrons are in 
thermodynamic equilibrium with a thermostat. Therefore, electron distribution under the condition (El) can be 
described by the conventional Fermi-Dirac function, where the energy of “bare” electron should be replaced with the 
energy of dressed electron da,)©. Substituting both this Fermi-Dirac function and the scattering probability © 
into the conventional kinetic Boltzmann equation, we can analyze the stationary (dc) transport properties of dressed 
electrons in graphene. Within this approach, we take into account the two key physical factors arisen from a dressing 
field: (i) modification of the electron energy spectra m and © by the dressing field; (ii) renormalization of the 
electron scattering probability ©-(nni) by the dressing field. 


RESULTS AND DISCUSSION 


Let us focus our attention on the dc conductivity of the dressed electrons. Generally, the density of the conduction 
electrons can be tuned by applying a bias voltage which fixes the Fermi energy, £^7-, of electron gas 0. Assuming the 
Fermi energy to be in the conduction band and the temperature to be zero, let us apply a stationary (dc) electric field 
E = {Ex,Ey) to the graphene sheet. It follows from the conventional Boltzmann equation for conduction electrons 
(see, e.g.. Refs. lilHl) that the electric current density, J, is given by the expression 

J = ^ y* [E • v(k)] r(k)v(k)^(£:k - eF)d^K (12) 

k 


where v(k) = (l//i)Vk£k is the electron velocity, and r(k) is the relaxation time. In the most general case of 
anisotropic electron scattering, this relaxation time is given by the equation 


1 

T(k) 


E 1 


T(k')E • v(k') 
r(k)E • v(k) 


Wk'k- 


(13) 


Substituting the scattering probability of dressed electron © into Eq. ED, we can obtain from Eqs. E1)-(IIS1) the 
conductivity of dressed graphene, = Ji/Ej. 

To simplify calculations, let us consider the electron scattering within the s-wave approximation [^ . where the 
matrix elements C/k'k do not depend on the angle 0k'k = (k',k)- Substituting the probability © into Eqs. (IT^ - (IT^ . 
we arrive at the isotropic conductivity of a graphene dressed by a circularly polarized field, ctc = ctxx = cTyy, which is 
given by the expression 


^ ^ 1 - £g/4£f 
CTo 1 + 3£2/4£| ’ 


(14) 


where ep >£g- It is seen in Fig.[^ that the conductivity (IT4l) monotonically decreases with increasing field intensity 
/o = eo£'gc/2. Physically, this behavior is a consequence of decreasing Fermi velocity, Vi7(k) = (l//i)Vk£k|e=£F j 
with increasing field amplitude Eq (see Fig. [5^). For the case of a dressing field linearly polarized along the x axis, 
the conductivity is plotted in Figs. ©i-©;. There are the two main features of the conductivity as a function of the 
dressing field intensity: Firstly, the conductivity oscillates, and, secondly, the giant anisotropy of the conductivity, 
<yxxl<^yy ~ 10 appears (see the insert in Fig.©). The oscillating behavior arises from the Bessel functions which 
take place in both the energy spectrum © and the scattering probability ©. As to the conductivity anisotropy, it 
is caused by the field-induced anisotropy of the energy spectrum ©. Namely, the linearly polarized dressing field 
turns the round (isotropic) Fermi line of unperturbed graphene into the strongly anisotropic ellipse line (see Fig.©). 
As a result, the Fermi velocities of dressed electrons along the x,y axes are strongly different and the discussed 
anisotropic conductivity appears. It should be stressed that the aforementioned features of electronic properties are 
typical exclusively for linear electron dispersion and, correspondingly, do not take place in a dressed electron gas 
with parabolic dispersion [s^. To avoid misunderstandings, it should be noted also that the zeros of conductivity in 
Fig. © lie within physically irrelevant areas pictured by dashed lines. Formally, these irrelevant areas correspond to 
the broken condition w ^ vkp^ which is crucial for the correctness of the energy spectrum © at the Fermi energy. 
Thus, the zeros have no physical meaning and should be ignored. 
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FIG. 3: The conductivity of dressed electrons in graphene for the dressing field with the different polarizations: (a) circularly 
polarized dressing field; (b)-(c) dressing field polarized along the x axis. Physically relevant regions of the field parameters, 
where the developed theory is applicable, correspond to the solid lines in the plot (b) and wide areas between the dashed lines 
in the plot (c). 


It is seen in Fig. [3] that the behavior of conductivity is qualitatively different for the dressing field with different 
polarizations. Physically, the strong polarization dependence of electronic transport follows directly from the strong 
polarization dependence of energy spectrum of dressed electrons. Namely, the energy spectrum of electrons dressed 
by a circularly polarized field ([2]) is isotropic and has the field-induced gap ([3|) at the Dirac point. In contrast, the 
energy spectrum of electrons dressed by the linearly polarized field ([5]) is gapless and has the field-induced anisotropy 
arisen from the Bessel function in Eq. O- These differences in the spectra m and dni) lead to the discussed difference 
of transport for electrons dressed by circularly polarized light and linearly polarized one. It should be noted that 
an electromagnetic field can open energy gaps within conduction and valence bands at electron wave vectors k 0 
(see, e.g.. Refs, [s^ - ldlj ]). These gaps arise from the optical (ac) Stark effect and take place at resonant points of the 
Brillouin zone, where the condition of w = 2vk is satisfied. Certainly, the basic expressions (l2|)-([7|) are not applicable 
near the Stark gaps. However, these gaps lie far from the Dirac point in the case of high-frequency dressing field. 
Therefore, they do not influence on low-energy electronic transport under consideration. 


CONCLUSIONS 

We have shown that the transport properties of electrons in graphene are strongly affected by a dressing field. 
Namely, a circularly polarized dressing field monotonically decreases the isotropic conductivity of graphene, whereas 
a linearly polarized dressing field results in the oscillating behavior of the conductivity and its giant anisotropy. As a 
result, the dc transport properties of graphene can be effectively controlled by a strong high-frequency electromagnetic 
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field. From the viewpoint of possible applications, the discussed effect can make graphene more tunable. Particularly, 
the switching times for conductivity of graphene controlled by a high-frequency held are expected to be shorter then for 
the case of conventional electrostatic control of conductivity by gate electrodes. This can create physical prerequisites 
for novel graphene-based optoelectronic devices. 
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SUPPLEMENTARY INEORMATION: 

Full and consistent derivation of the energy spectrum and wave functions of dressed electrons in graphene 


Let us consider a graphene sheet which lies in the plane {x, y) at z = 0 and is subjected to an electromagnetic 
wave propagating along the z axis (dressing electromagnetic field). Then electronic properties of the graphene are 
described by the Hamiltonian 


a = vcr ■ (/ik — eA), (15) 

where cr = ay) is the Pauli matrix vector, k = ky) is the electron wave vector in the graphene plane, v is the 
electron velocity in graphene near the Dirac point, e is the electron charge, and A = {A^^Ay) is the vector potential 
of the electromagnetic wave in the graphene plane. Solving the Schrodinger equation with the Hamiltonian (I15L we 
can find the energy spectrum of dressed electrons and their wave functions as follows. 


I. Circularly polarized dressing field 


For the case of circularly polarized electromagnetic wave, its vector potential A = {Ax, Ay) can be written as 


A = I — cos U!t, — smwr 
w a; 


(16) 


where Eq is the electric field amplitude of the wave, and w is the wave frequency. Then the Hamiltonian (1151) is 

n^iio + nk, (17) 

where 

vcEq 


no={- 


UJ 


0 

jiujt Q 


and 


Uk = 


0 

v{hkx + ihky) 


v{hkx — ihky) 


0 


The nonstationary Schrodinger equation with the Hamiltonian (1181) . 


(18) 


(19) 


( 20 ) 


describes the time evolution of electron states in the Dirac point (k = 0). The exact solutions of the equation (l20l) 
can be sought in the form 


^0 = 


( 21 ) 


where $2 2 (r) are the basic functions of the 2x2 matrix Hamiltonian (I15p . and a, A and B are the sought parameters. 
Substituting the wave function (^11) into the Schrodinger equation (1^ . we arrive at the system of algebraic equations 


A 

A 


fUjJ 


( veEo\ 


\ UJ J 

The condition of nontrivial solution of the system 

a + 




= 0, 


0 

0 . 


( 22 ) 


huj veEo 


veEp _ hoj 

UJ ^ 2 
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gives the two different parameters, a = ±0/2, where 


^ (23) 

Therefore, there are two sets of solutions of the system (IdlT) . which correspond to these two parameters and satisfy 
the normalization condition, \A\'^ + |i?p = 1. As a result, there are two wave functions (I^Tl) . 




O ± feu 
20 




fl —|— feiU 
20 


$'2(r)n 


iujtj2 


(24) 


which exactly describe electron states of irradiated graphene in the Dirac point (k = 0). Since the two wave functions 
(12411 are the complete function system for any time t, we can seek the solution of the Schrodinger equation with the 
full Hamiltonian (1171) as an expansion 


V'k = a^{t)'ip^ + a (<)V’o ■ 

Substituting the expansion (1^51) into the Schrodinger equation with the full Hamiltonian (1171) , 

ot 

we arrive at the system of two differential equations for the coefficients a^(t), 

e 


ia/it) = 


1^ 

O 


{kx COBujt + kyS\nujt)d^[t) — [kx —iky) ^® 


± [kx ± iky) 


n — hu 
20 




it) 


ia (t) = 


Wt 

”o 


^ {kx cos ujt + kySmojt)a (t) + {kx + iky) ^ ‘^'^*a'^{t) 


- {kx - iky) 


O — feu 
20 


j{n/h+uj)t + 


it) 


(25) 


(26) 


(27) 


where 


2veEo 
Wo = - 

U! 

is the characteristic kinetic energy of rotational electron motion induced by the circularly polarized field. In what 
follows, we will assume that the field frequency cu is high enough to satisfy the condition 


Wo / feu -C 1. 


(28) 


Then we have 

Wq ^ Wo ^ 0±fecu^ fl — huj 
”0” ^ fe)7 ^ ’ 20 ^ ’ 20 

and, correspondingly, Eqs. (I?7l) take the form 


(29) 


ia+{t) = -v-^{kx-iky)e-^^^-^‘^'>*/'^a-{t), 

id-{t) = -v^^{kx+iky)E^^-'^^'>^/^a+{t). (30) 

Let us seek solutions of Eqs. o as af{t) = C^e *(±t2/2T^i‘^/2+ek)t/fi^ where and are the sought parameters. 
Solving the system of equations O, we arrive at 


£k = (£ff/2)^ ± {hvk)"^, 


(31) 
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and 


Eg = W {hw) + ( 


f 2veEQ \' 

w J 


— huj. 


(32) 


Following the conventional terminology, Eq. ()31|) describes the spectrum of quasienergies of electrons in graphene 
dressed by a circularly polarized field, where signs “+” and ” correspond to conduction band and valence band, 
respectively. The two wave functions (12511 . which correspond to the two energy branches m, can be written as 


= e 


_ „-iSkt/h 


±-1 


kkl T ^g/2 


2kk| 


g/^^-ie/2 


e ft + hw . 




2n 


-$'(r)e" 


2|ek| IV 212 |e|V 212 ^ 


(r) 


(33) 


where 9 is the azimuth angle of the wave vector, k = {kcos9, fesind). Let us write the basis functions 2 ( 1 ') 
conventional Bloch form, $' 12 ( 1 ") = <l’i^ 2 (r)i^k(r), where $i 2 (r) are the periodical functions arisen from atomic tt- 
orbitals of the two crystal sublattices of graphene, </5k(r) = /\fS is the plane electron wave, and S is the graphene 
area. Keeping in mind inequalities (I29L we arrive from (1331) to the sought wave functions of dressed electrons in the 
final form. 


^k = ^k(r)e— 


(34) 


II. Linearly polarized dressing field 

For the case of electromagnetic wave linearly polarized along the x axis, its vector potential A = {A^,Ay) can be 
written as 


A = ( — coswf, 0 1 . 

OJ 


(35) 


Then the Hamiltonian (fT51) is 


'H='Ho+nk: 


(36) 


where 


no = 


0 —1\ vcEq 
-I 0 ~ 


■ cos tot. 


(37) 


and 


Uk = 


0 v{hk^ — ihky) 

v{hkx + ihky) 0 


(38) 


The nonstationary Schrddinger equation with the Hamiltonian (1371) . 

th^- = HoVo, 
at 


(39) 


describes the time evolution of electron states in the Dirac point (k = 0). The exact solutions of this Schrddinger 
equation can be sought in the form 


ipo = [-4$((r) + B4>2(r)] exp 


. veEo . 

—la ■ ^ smcut 


(40) 
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where a, A and B are the sought parameters. Substituting the wave function (|40)) into the Schrodinger equation (l3^ 
with the Hamiltonian (I37L we arrive at the system of algebraic equations 


Aa + 5 = 0 
A + Ba = 0. 


The condition of nontrivial solution of the system dUD, 


a 1 

1 a 


= 0 , 


(41) 


gives the two different parameters, a = ±1. Therefore, there are two sets of solutions of the system dm, which 
correspond to these two parameters and satisfy the normalization condition, + |5p = 1: 


.4 = -B=T ,42, 

As a result, there are two wave functions (HOI) . 


'<Pt 


[$;(r) ± $^(r)]exp 


,.veEo . 

±l , ^ SlULUt 


(43) 


which exactly describe electron states of irradiated graphene in the Dirac point (k = 0). Since the two wave functions 
(I43p are the complete function system for any time t, we can seek the solution of the nonstationary Schrodinger 
equation with the full Hamiltonian (1361) as an expansion 

-0k = a+(t)-0,)' + (44) 


Substituting the expansion (H41) into this Schrodinger equation with the full Hamiltonian (1001) . 

(45) 

ot 

we arrive at the system of two differential equations for the coefficients a^(t). 




(46) 


It follows from the Floquet’s theorem that the functions a^{t) can be written as 

a±(t) = e-*''‘‘*/^a±(t), (47) 

where Sk is the electron quasienergy in the irradiated graphene (the energy spectrum of dressed electrons), and a^(t) 
are the periodical functions with the period T = 271 juj. Let us apply the Jacobi-Anger expansion, 

C30 

n— — oo 

to the exponents in the right side of Eqs. (HOI) and expand the periodical functions a^(t) into the Fourier series 

OO 

S±(t)= ^ c±e—0 


n— — oo 
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Then the differential equations (l46l) can be transformed to the algebraic equations 

2veEo \ 
huj^ ) 


^ vkx -^ — nuj 

\ h 


+ T 


OO / 

) c+ + ivky f 

OO 

+ nuj] c“ + ivky Y 


= 0 , 


^n-m^rn 




(48) 


where Jm{z) is the Bessel function of the first kind. Equations (H51) can be easily solved in the case of high-frequency 
field satisfying the conditions 


kiv ^ Ek, uj ^ vk. 

In this simplest case, Eqs. (H51) for n ^ 0 can be reduced to the equation 


c 


± 

n 



m— — oo 


( 2veEo \ 


(49) 


(50) 


Keeping in mind that |c^| < 1 and \ Jn[z)\ < 1, Eq. (1501) leads to the evident solution ~ 0. After substitution of 
this solution into Eqs. (l48l) . the expressions (l48l) are reduced to the two equations 

{yk^ - + ivkyC^ Jo = 0 ’ 

(vk^ + y)^o + Jo = 0 - 

(51) 

Solving the system of the two algebraic equations (ISTI) accurately, we can easily obtain both the coefficients and 
the energy spectrum of dressed electrons, 


£k = ±hvkf{e), 


(52) 


where 


/(6») = y cos2 9 + J§ sin^ 

Correspondingly, the sought wave functions of dressed electrons (|44)) take the final form 
V’k = 

sin 6 ( 2veEo \ 


[$i (r) ± ^r)] 

(^) [<hi(r)T‘i>2(r)]e=F*(-^“/'‘-^)--A 


(53) 


Taking into account Eq. (l52l) . the two conditions (l4^ can be reduced to the solely condition. 


u! ^ vk 

which describes borders of applicability of Eqs. (15^ and (1551) . 


(54) 













